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Abstract 


One of the major applications of the Domain Decomposition Time March- 
ing Algorithm is the coupling of the Navier-Stokes systems with Boltzmann 
equations in order to compute transitional flows. Another important ap- 
plication, is the coupling of a global Navier-Stokes problem with a local 
one in order to use different modelizations and/or discretizations. Both of 
these applications involve a global Navier-Stokes systems with non standard 
boundary conditions. The purpose of this work is to prove, using the classical 
Leray-Schauder theory, that these boundary conditions are admissible and 
lead to a well posed problem. 
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1 Introduction 


In this paper we study the Navier-Stokes equations with the new boundary 
conditions introduced by the application of the Domain Decomposition Time 
Marching Algorithm ( [6], [9], and [10]). These boundary conditions are of 
slip type , they could appear either by the application of the Domain Decom- 
position Time Marching Algorithm to a Navier-Stokes problem ([9] and [10]), 
or to a Navier-Stokes/Boltzmann coupling ([7] and [9]). In the latter case, 
they are similar to the analytical slip boundary conditions introduced in [5], 
which are derived from kinetic theory in order to replace solving Boltzmann 
equations by solving Navier-Stokes equations in the transitional regime. 

We study here only the stationary problem. The treatment of the time 
dependent problem will follow using the same ideas developed here (see [9]), 
and the classical proofs for the standard boundary conditions (see [l]-[3], 
and [4]). We begin by describing the strong formulation of the problem in 
the first paragraph, then we set the preliminary results necessary for our 
study. In the third paragraph we show the equivalence of the strong and the 
weak formulations under some regularity hypothesis, from which we deduce 
the admissible boundary conditions. We present, then, the study of the 
stationary problem. Paragraph 6, deals with the uniqueness results when 
the data are sufficiently small. Finally, we present, in the last paragraph, 
some conclusions. 


2 Motivations 

2.1 The general coupling strategy 

For coupling external Navier-Stokes equations, with local Navier-Stokes equa- 
tions (dense regimes) or local Boltzmann equations (transitional regimes), we 
introduce two domains, a global one 0, a local one Qy included in 0, and 
an interface T; (Fig. 1 in which T e denotes Too). The global solution W on 
0 and the local solution Ui oc on Qy, are matched by the following boundary 


1 



Figure 1: The global geometry 


conditions, inspired of Schwarz overlapping techniques : 

W = given imposed value on T (;o , 

n ■ a(W) ■ t = n ■ cr(Ui 0C ) ■ r on the body T 0 , (equality of friction forces) 

< q(W) ■ n + n ■ a(W) ■ v = q(Ui 0C ) • n on T 0 , (equality of total heat fluxes) 
v ■ n = 0 on r 0 , 

Uloc = 0 on r 0 , Uioc = W on the interface T;. 

Above, n ■ a ■ n and n ■ a ■ r respectively denote the normal and the 
tangential force exerted by the body on the flow, with n the unit normal 
vector to the body oriented towards its interior. 

The calculation of Ui oc and W satisfying the above boundary conditions 
is then obtained by the Domain Decomposition Time Marching Algorithm, 
which was introduced by Le Tallec and Tidriri ([6], [9] and [10]) and which 
leads to the following algorithm : 

Initialization 

1. Guess an initial distribution of the conservative variable W in the 
global domain 0 ; 

2. Advance in time this distribution by using the global Navier- Stokes 
solver on Ni time steps, with Diri chief type boundary conditions on the body 

r 0 ; 
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3. Deduce from this result an initial distribution of the local variable Ui oc 
on the interface I\- and in the local domain fly ; 

4. Advance in time this distribution by using the local solver on N 2 time 
steps with Dirichlet boundary conditions on I\- and T 0 . 

Iterations 

5. From Ui 0C} compute the friction forces n ■ cr(Ui 0C ) • r and heat flux 
q(Ui 0C ) ■ n on the body T 0 ; 

6. Advance the global solution in time (iVi steps) by using the global 
Navier- Stokes solver with the above viscous forces as boundary conditions 
on r 0 ; 

7. From W, compute the value of Ui oc on the interface T; ; 

8. Using this new value as Dirichlet boundary conditions on Ty advance 
the local solution in time (iV 2 steps) and go back to step 5 until convergence 
is reached. 

This algorithm completely uncouples the local and the global problems 
which can therefore be solved by independent solvers. 

A parallel version is also quite possible although it is generally wiser to 
use parallel solvers within steps 6 and 8. 

Remark 2.1 The local problem can be either Navier-Stokes or Boltzmann 
equations (see [6]- [10]). 


2.2 The global Navier-Stokes problem 

The global domain 0 is discretized using node centered cells defined on an 
unstructured grid. Then, at each time step n and for each cell z, we solve 


W n+1 - W r ‘ 


ic, 


At 


+ £ 


.^ T .. ^ 8C t ndC, 

jev(t) j 


F C (W 


71 + 1 ^ 


• Til 


+ F D (W n+1 ) ■ rii -\- / F(W n+1 ) ■ rii = — / F 0 -n % . 

Jdc.-T Jdc.nr^ Jdc t nr 0 

The fluxes Fc and are computed at time step n + 1 and linearized, with 
for example Fc computed by an Osher approximate Riemann solver [9], and 
[ 10 ]- 
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Figure 2: A boundary cell 

On the body T 0 , because of our special choice of boundary conditions, 
the flux is given by 

° \ 

Hi ■ a(W n+1 ) ■ n i 

n-i ■ o-{Ui oc ) ■ n 
V -q{Ui 0 C )-n t ) 

where the aspect of a boundary cell C; is described in Fig. 2. 

In other words, friction forces and heat flux are given explicitly as pre- 
dicted by the local solver and the mass flux is imposed to zero. Then, in 
order to have a well-posed problem, at least in the incompressible case (see 
next paragraphs), the normal stress (the multiplier of the zero mass flux con- 
straint) cannot be imposed and must be obtained from the solution W n+1 . 



Remark 2.2 Imposing friction forces to the global solution instead of no 
slip boundary conditions allows to have an accurate solution away from the 
boundary layer even with a coarse mesh (see [9]). 


3 Strong formulation 


Let 0 be the domain occupied by the fluid, T its boundary as described in 
figure (1); we assume that 0 satisfy : 

0 is a simply connected bounded domain in IR^, 

1 = dfl is of class if 1 , 


( 1 ) 


r 6 U = r and h n I# = I, 

T& and T (>CJ are compacts of nonzero measure. 
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In the steady incompressible case, the global problem in the coupling 
strategy described above consists in finding u and p satisfying the following 
equations : 


z/A u + u'Vu + Vp = / 

in 

0, 

(2) 

divu = 0 

in 

0, 

(3) 

u ■ n = 0 

on 

r 6 , 

(4) 

u = h 

on 

Too, 

(5) 

z/(Vn • n)r = g T 

on 

r 6 , 

(6) 


where n is the unit normal to T, r is any unit normal to n. 

We look for weak solutions of problem (2)-(6), i.e. (2), (3) are satisfied 
in a distribution sense and equations (4)- (6) are satisfied in some Sobolev 
spaces. 

Remark 3.1 From the remark (2.1) we observe that g T can be issued from 
a local kinetic model (Boltzmann equations) or a local continuous model 
(Navier-Stokes system). 


4 Preliminary results 

Let H m ~ 2 (T) be a space of trace functions in fL m (0) with the following 
norm : 

IHIm-iy = inf{||u|| mi n/u £ H m (kl),u\ r = f}, (7) 

and V the closed subspace of (77 1 (0)) JV defined by 


V = {u £ (H 1 (Ct)) N \divv = 0, v ■ n = 0 on r&, v = 0 on Too}, (8) 
and define 

IMI = ||Vu||o, 2 ,fi. (9) 

We have the following classical lemma. 
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Lemma 4.1 ||.|| is a norm ofV equivalent to the norm ||.||i, 2 ,fi- 
Lemma 4.2 ([11]) Under hypothesis (1), we have 


n : T — 

* IR V 

(10) 

x — 

■> n(x) 


is C 1 and hence, the trace operator 



- 

-> Hi(T) 

(11) 

u — > u 

■ n 


is continuous. 



Lemma 4.3 ([14]). Let 0 satisfy (1). We have : 
The injection 

W hp (fl) C L qi ( 0) 


is compact for qi satisfying : 



1 < qi < oo 

if p > N 

(12) 

1 < qi < q with q given by — - 

q 

= -- 4 */ 1 <P< N. 

p iv 

(13) 

In particular we have for p = 2 : 



H^fl) C L q (fl) for 

2 N 

2 < q < 

~ H ~ iV-2’ 



the injection being continuous. For N = 2, this injection takes place for any 
finite q such that q > 2. 


Lemma 4.4 ([12]) Let ft be a bounded connected domain. Then 
(i) The operator grad is an isomorphism from Lq(TI) into V° , 
(u) the operator div is an isomorphism from (V 0 ) 1 into Lq(Q) 
where Lg(fi) and V° are defined by 

v° = {he H~ 1 (n) N ,<h,v>= o, Vu g v}. 

L 2 Jtt) = {he L 2 Jtt), [ h(x)dx = 0}. 
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We now set 


N 


dv ; 


b(u,v,w) = y / Ui— — -w;dx,u,v,w £ V 

‘dx, 


(14) 


< /, v >= f fvdx,v£V (15) 

Ju 

<g T ,v> r b = / g T r ■ vd^,v eV (16) 

Jr b 

((u,n)) = J 'Vu'Vv dx } u } v £ V (17) 

(u,v) = J uvdx 7 m,d £ K (18) 

If / is in (i7 _1 (0)) JV and g T in all the above operators (14)-(18) 

are well posed. Moreover, the operator defined by (16) and the trilinear form 
defined by (14) satisfy the following properties : 


Lemma 4.5 The mapping from V into IR defined by (16) is continuous. 


Proof: The proof follows directly from lemma 4.2. 

Lemma 4.6 (i) The trilinear form b is continuous on V 2 X (V fl (L N (Q)) N 


(u) Tet u,v,w G V . Then we have the following relations : 

b(u } v } w) = —b(u } w } v) } (19) 

b(u , v, v) = 0. (20) 

(in) Tet {u n } n > 0 be a sequence in V such that u n u weakly in V . Then : 

lim b(u n} u n . v) = b(u } u. v), Vv £ V. (21) 

71^00 V ' V ' V ' 


Proof 

We present only the proof of (i) and (ii). For the proof of (iii), see [14] or 


[13]- 


(i) Let u G V then U{ £ H 1 (fl) ; from lemma 4.3, we deduce : 


Ui £ T 9 (Q), 


1 

6 


1 1 
2 ~ N’ 
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which implies 


Ov ■ Ov ' 

^ 1 “ W^dx | ^ 1 1 1 1 0,q,Q 1 1 Tj 1 1 0,q,Q \ \ \ \ 0,N,Q 


'Jn J dx 3 1 1 - " 

— ^IkjlklNklMk- 

Therefore we conclude to our result. 

(ii) By dehnition 

b(u,v,w) = / UiVjiW 

Ju 

After integration by parts, we obtain 


( 22 ) 




b(u } v, w) = - Vj(ui w j) ti + J 


u ■ nv ■ w 


= — VjUiWji — / VjWjUi i + / u ■ nv ■ w. 

Jfi Jfi ’ JdQ 


Since we have it -re = 0 on dfl and ipy = divu = 0 on 0. Therefore, 

we get 

b(u,v,w) = — VjUiWji = —b(u,w,v). 

Ju 

Remark 4.1 For N < 4 ; it follows from lemma f.3 that V fl (L N (0,)) N = V 


5 Existence Result 

We now go back to our initial formulation: 
Find u and p such that : 


is A u + u ■ Vu + Vp = / 

in 

0, 

(23) 

divu = 0 

in 

0, 

(24) 

u • n = <f> 

on 

Ik 

(25) 

u = h 

on 

Too, 

(26) 

isfVu • n)r = g T 

on 

IV 

(27) 
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5.1 Hypothesis on cf) and h 

The data </> on r& and h on Too must verify some regularity and compatibility 
properties. More precisely we assume that : 



(28) 


(29) 

f h ■ n = 0, 

(30) 

Jr co 

-e- 

II 

o 

(31) 

J r 6 

Lemma 5.1 The relations (28)-(31) imply the existence 
such that : 

of u € 

V.u = 0 in 0, 

(32) 

u ■ n = <f> on r & , 

(33) 

u = h on Too. 

(34) 

Proof 

Let if; be a function defined on T such that : 


( cf> ■ n on r & , 



(35) 

{ h on Too. 



From the assumed hypothesis, on T, </> and h we get : 


Therefore, from the trace theorem, we can take u £ H x (Vl) N such that : 
rt|r = V’- By using the relations (30) and (31) and the Green Formula, we 
obtain 

/ V • u = 0. 

Jn 
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Using the regularity of it, we have V • it £ Uq(Q). 

Lemma 4.4 shows that there exists v £ H[( 0) such that V • u = V • it. 
Finally, by taking now u = u — v, we can see easily that u verify the relations 
( 3 2 )-( 34 ). 

Now, we are going to pick out some function ( in H 1 ( 0) such that : 

u = rot (. ( 3 6) 

To achieve this, we will use the following proposition : 

Proposition 5.1 ([13]) Under the regularity hypothesis (1) and for N = 3 
ice have 

rotH 1 ^ 0) = {it £ L 2 (Q) N , divu = 0, f u ■ n = f u ■ n = 0} ( 3 7) 

Remark 5.1 This result is also valid for N = 2 ([13]). 

Under the hypothesis ( 3 0)-( 3 l) and the relations ( 3 2)-( 3 4) we have : 



Consequently, from Proposition 5.1, there exists ( £ iL 1 (0), such that we 
have ( 3 6). From now on we assume that N < 3. 


5.2 Weak formulation 


Let u £ V. After multiplication of (2 3 ) by u, integration by parts and by 
taking into account the boundary conditions (25)- (27) and the relation (24) 
we obtain the following variational formulation : 


( Find it 0 £ V such that 


< 


v((u 0 , v)) + b(u , it 0 , v) + b(u 0 , u, v) + b(u 0 , it 0 , v) 


= — n((h, u)) — b(u } u, v) + 
Vn £ U C 7^(0)^ 


f 9 tT ■ v + [ fv, 
r b 


(40) 
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5.3 Equivalence between the two formulations 

We will show here that the variational formulation (40) is equivalent to the 
strong formulation. This is the goal of the following proposition : 

Proposition 5.2 If (u,p) is a smooth solution of (2)-(6), then it satisfies 
(40) ■ Conversely, if u 0 is solution of (40), then there exists a unique function 
(up to a constant) p £ L 2 (fl) such that u 0 + u and p satisfy : 

(i) (2) and (3) in a distribution sense, 

(ii) (4) in H*(T h ), 

(iii) (5) in 

(iv) (6) in 
Proof 

The direct theorem was showed in our introduction of the weak formula- 
tion. We then show only the converse. 

Let 

V = {f £ (H((tt)) N / dim) = 0}. 

Let u 0 be a solution of problem (40) then after multiplication by a test 
function in and integration by parts and density, we get 

< — rAii + (it • V)u — /, f > = 0, \/f) £ V (41) 

where <,> denotes now the duality between the spaces ’D'(Jl) N and 
V(fl) N . We consider the operator —grad{= —V) £ H -1 (fl) N ), 

and IZ(-grad) its image in iL -1 ( fl) N . We have the following lemma ( see 

[ 12 ]) 

Lemma 5.2 7 Z(-grad) is identical to V° with V° defined by : 

v° = {he < h,f>= o,w> £ v}. 
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The application of this lemma implies the existence of p £ Xq(Q) such 
that 

-vAu + u- Vu-f + Vp = 0 in 
As it £ V, it • Vu and / £ X 2 , then 


— vAu — 'Vp = div( — v'Vu + pld) 

= (f-u-Vu)e(L 2 ( n)) N . 

Applying the Green formula implies that there exists 
q = ( — vVu + pld) ■ n £ 

such that 

/ — div( — rVti + p/d) • u 

Jo. 

+ / j/VtiVt) — / pdivv = q ■ v, Vu £ (AT 1 (0) ) 

By replacing div( — v'Vu + pld ) by (/ — u ■ Vrt) we get 


JV 


n((n, n)) + b(u } it, v 


/ pdivv 
Ju 



Vn £ (^(fi))^ 


( 42 ) 


By subtracting (43) from (40), it remains 



i i 

Now let w arbitrary element of (T 6 ) , where -f/oo(r&) is defined by : 


tf 0 2 0 (T 6 ) = {n £ L 2 (Tb); 3 w £ i7 1 (0), w = 0 on Too andic = v on Tf, } 
(see [11] for more details about this space). 
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Let v £ V such that 


v ■ n = 0 on r & , 
v ■ t = w on r & , 
v = 0 on Too. 

Applying the precedent equality, it remains 


/ q ■ tw = / g T w,Vw £ H£ 0 (T h ), 
>r b Jr b 


therefore 


q ■ t = t ■ (-rVu + pld ) • n 
= g T in H~^(T b ). 

Remark 5.2 We cannot impose the normal stress in (3)-( 7) because of the 
presence of a pressure term in q. 

5.4 Existence result 

Now, we are able to set out the main result of this section. 

Theorem 5.1 Let f £ ; then there exists u 0 £ V solution of (40). 

Proof 

We proceed as in the usual proof of existence for the incompressible 
Navier- Stokes equations. Let us define 

[if m (u)i u] = z/||u|| 2 + n((u, it)) + &(u, ft, u) — L g T ru 


— fu + b(u , ft, 

Ju 


By using the continuity of b we obtain : 
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( 43 ) 


with q satisfying (12)-(13). In addition we have : 

b(u, u, u) < c 2 ||'u || 2 2 n ||rt||. (44) 

Using lemma 4.1 and lemma 4.2 we arrive to : 

[ip m (u),u\ > n||n|| 2 + b(u, u, u) 

“IMIMMI + C 2|| M lll,2,fi + c 3 ||^T || _ + ||/||-l,2,fi)- 

So, to be able to use the Brouwer fixed point theorem, it is sufficient to 
show that we have for some f3 > 0, the following relation : 

n||n || 2 + b(u } u } u) > /3||m|| 2 ,Vm £ V. (45) 

This is the goal of the following lemma : 

Lemma 5.3 V 7 > 0 we can choose u verifying (32)-(33) such that 

|&(rt, u, u)\ < 7 ||n|| 2 . (46) 

From the above lemma, we conclude to the existence of a given k such 
that 

[fm{u), (u)] > 0, 

for any u £ V with ||rt|| = k. Then the Brouwer hxed point theorem can be 
applied to the function t/j m (u) inside any ball of radius k belonging to any 
given finite dimensional approximation V m of V . From this, we deduce the 
existence of a Galerkin approximation u m £ V m of the solution u of (40). By 
standard compactness arguments, it follows that u m weakly converges to a 
solution u 0 of (40) in V . (see [5] for more details). 
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6 Uniqueness results 

Under the assumption that N < 3, we have the following result : 
Theorem 6.1 Assume that u in (L N (Q)) N is sufficiently small so that 

\b(v } u } v)\ < pn|| 2 ,Vn £ V, 
and v is sufficiently large so that 


v 2 > 4C(\\f\\v + llfl'rll-i.rj 
where C is the constant in (22) and 

d 

f = z/A u — E, — — u + f 

3 3 dxj 

and 


9t 


d „ 

n—u + g 
On 


T 


Proof 

We proceed as in the proof for the standard homogeneous boundary 
ditions (see for instance [13]). 

Suppose that u 0 is solution of the problem (40). Taking v = u 0 in 
and using (20), we obtain : 


V\\U 0 \ y — 


-b(u 0 , u, u 0 ) - v((u, u 0 )) - b(u, u, u 0 ) + / g T r-u 0 + fu 0 , 

•j r h j 


Using the Green formula, we obtain 


d 

Hkllt = — 0, ,, o)+ < nAO — S, Uj- U + f,U„ > 

Ox ; 


d 

+ < v-TT-u +9 tO“o >r t 


= -b(u 0 ,u,u 0 )9- < f,u 0 > + < 9rOX. 0 >r b 


(47) 

(48) 

con- 

140) 

(49) 

(50) 

(51) 

(52) 
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Since / G (H-^nyf and u G (Hff Sl)) N , A u G {H-fltiff and then 
/ G (i? _1 (0)) Jv . Similarly, since g T G (-fT - ^ (T^ ) ) iV and u G (i? 1 (0)) JV we 

have — G (H~^ (TAl^ and then o T G (H~^ (Tb)) N , and the equality above 
on 

is well dehned. The rest of the proof is an adaptation of the proof for the 
standard homogeneous boundary conditions (see for instance [13]). 

7 Conclusion 

In this paper we have shown that the slip boundary conditions resulting 
from the application of the Domain Decomposition Time Marching Algo- 
rithm to either Navier-Stokes/Navier- Stokes coupling ([9] and [10]) or Navier- 
Stokes/Boltzmann coupling ([9], [7], and [8]) are admissible and lead to a well 
posed global problem at least in the incompressible case. 
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